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Abstract 



We search for three-neutron resonances which were predicted from pion dou- 
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ble charge exchange experiments on 3 He. All partial waves up to J = | are 



nonresonant except the J 71 " = | one, where we find a state at E = 14 MeV 
energy with 13 MeV width. The parameters of the mirror state in the three- 
proton system are E = 15 MeV and T = 14 MeV. The possible existence of an 
excited state in the triton, which was predicted from a H( 6 He,a) experiment, 
is also discussed. 

PACS numbers: 21.45. +v, 24.30.Gd, 27.10. +h, 25.80.Gn 
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I. INTRODUCTION 



One of the clearest indication of the spin dependence of the nucleon-nucleon interaction 
is the fact that there exists a bound triplet deuteron, but not a singlet deuteron or dineutron. 
The possibility that by adding one or two more neutrons to an unbound dineutron one might 
wind up with a bound tri-neutron or tetra-neutron has long been discussed. However, precise 
experiments and calculations show no indication for the existence of such bound structures 
|Ij]. Nevertheless, the existence of resonant states of these nuclei, which have observable 
effects, cannot be excluded. 

Sperinde et al. found that the differential cross section of the 3 He(7r~, TT + )3n pion 
double charge exchange reaction strongly differs from what one would expect from a pure 
phase-space description of the final state. The discrepancy seemed to be well explained by 
the assumption that there exists a three-neutron resonance at 2 MeV energy with 12 MeV 
width. The authors suggested that this state is probably a member of a T = | isospin 
quartet in the A = 3 nuclei, another member being the three-proton state found in || . The 
mechanism of the double charge exchange experiment was studied theoretically in |4j . The 
authors found that by taking into account the final state interaction between the outgoing 
neutrons, one can get good agreement with the experimental results of |fjj. However, the 
nature of this final state interaction, whether it produces a three-neutron resonance or not, 
was not studied in ||. 

Recently, a more thorough experimental study of the 3 He(7r~, n + )3n process was per- 
formed 0. The authors pointed out that the analysis of [0] was incorrect, because a four- 
body phase-space had been divided out of the cross section, instead of the correct three- 
body one. And this mistake leads to a peak in the cross section at low energies. However, 
three-neutron structures also appear in || at small angles. These structures at ~ 20 MeV 
three-neutron missing mass and with ~ 20 MeV width could be caused by three-neutron 
resonances. 

There are experimental indications of the possible existence of resonances in other A = 3 
nuclei, too. For example, in a recent H( 6 He,a) experiment the authors concluded that 
in order to explain their results, an excited state of the triton has to be assumed with 
E* = 7.0 ± 0.3 MeV excitation energy and T = 0.6 ± 0.3 MeV width. 

Theoretically the problem of a three-neutron resonance was studied by Glockle |7J , who 
calculated the pole position of the S matrix from the analytical continuation of the Faddeev 
kernel for three neutrons with pure 1 5'o N — N interaction. He found that the pattern of 
the pole trajectory ruled out the possibility of a low-energy | + three-neutron resonance in 
his model. 

The aim of this current work is to systematically investigate the possible existence of 
resonances in the lowest partial waves of the three-neutron system by using the complex 
scaling method. We shall also comment on the triton excited state. 

II. COMPLEX SCALING 

In coordinate space, resonance eigenfunctions, corresponding to the complex energy so- 
lutions of the Schrodinger equation 
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H\V) = (f + V)\V>) = E\V>), (1) 

show oscillatory behavior in the asymptotic region with exponentially growing amplitude, 
~ exp[i(« — ry)r] («, 7 > 0). Thus, the resonance eigenfunctions are not square-integrable. 
The complex scaling method (CSM) || reduces the description of resonant states to that 
of bound states, thus avoiding the problem of asymptotics. In the CSM we define a new 
Hamiltonian by 

H e = U(9)HU- 1 (9), (2) 

and solve the complex equation 

Ho\*o)=e\V e ). (3) 

Here the unbounded similarity transformation U(9) acts, in the coordinate space, on a 
function /(r) as 

U(9)f(r) = e^' 2 f{ve ie ). (4) 

(If 9 is real, U{9) means a rotation into the complex coordinate plane, if it is complex, 
it means a rotation and scaling.) In the case of a many-body Hamiltonian, (|4]) means 
that the transformation has to be performed in each Jacobi coordinate. For a broad class 
of potentials there is the following connection between the spectra of H and Hg ||: (i) 
the bound eigenstates of H are the eigenstates of Hg, for any value of 9 within < 9 < 
7r/2; (ii) the continuous spectrum of H will be rotated by an angle 29; (iii) the complex 
generalized eigenvalues of Hg, e res = E — iT/2, E, F > (where F is the full width at half 
maximum) belong to its proper spectrum, with square-integrable eigenfunctions, provided 
29 > I arg£ res |. These complex eigenvalues coincide with the S matrix pole positions. This 



method was tested for three-body resonances in |]10| , and was applied to the low-lying spectra 
of the 6 He, 6 Li, and 6 Be nuclei in |pT |. Further details and references of the method can be 
found there. 

Up to Eq. (H) our treatment of the three-body resonances is exact. The only approxima- 
tion we use here is that now we expand the wave function of Eq. ([3]) in terms of products of 
Gaussian functions with different widths. Thus, we select the square integrable solutions of 
Eq. @, which are the three-body resonances, and discretize the continuum. A term of this 
expansion looks like p\ exp[-(pi/7;) 2 ]r iimi (pi) • p 2 2 exp[-(p 2 /7 j ) 2 ]Y i2m2 (p 2 ), where h and l 2 
are the angular momenta in the two relative motions, respectively, and the widths 7 of the 
Gaussians are the parameters of the expansion. 



III. RESULTS 

We use the Minnesota effective nucleon-nucleon interaction \\L2\\ , together with the tensor 
force of fil||i~4|1 . The space exchange parameter of the tensor force is adjusted to the splitting 
of the 3 Pj phase shifts. As one can see in Fig. 1, this force produces p + p phase shifts which 
are in good agreement with experiment. The higher partial waves, not shown in Fig. 1, are 
also close to the (practically zero) experimental values. In the three-neutron wave functions 
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all n + n partial waves up to / = 2 are included. In some cases we checked the role of the 
/ = 3 partial waves and found them insignificant. 

Figure 2(a) shows a typical result of our calculations, in this case for J n = | . What 
one can see is the complex eigenvalue spectra of Eq. (|3|), with E = Re(e) and T = 2Im(e). 
The eigenvalues are all sitting on a half-life whose angle to the E axis is roughly 29. This 
line is the discretized rotated continuum of Eq. (^). If there were a resonance in this partial 
wave, it would be in the upper right triangle, separated from the continuum points. 

In Fig. 2(b) we show the result for J 71 " = | + , which is found to be the only resonant 
partial wave. The resonant state is shown by the circle. We find that the dominant terms 
in the resonant J n = | + wave function are the [(lil 2 )L, S] = [(11)1, |] and [(11)2, |] LS 
components, where L and S are the the total angular momentum and total spin, respectively. 
The other components have small but non-negligible effects. 

Because we solve Eq. (§) approximately, the resonance papameters depend slightly on 9. 
To find the optimal 9 value one should repeat the calculations with several different 9 values, 
and find a region where the resonance parameters are most insensitive to the change of 9. 



This way of finding the optimal 9 value is related to the complex virial theorem JT6|. To 
reduce computer time, we perform these calculations with a wave function which contains 
only the above mentioned two most important LS components. In Fig. 3 one can see that 
the optimal 9 angle is somewhere between 0.3 and 0.35 rad. Using such a 9 value in the full 
model the parameters of the resonance are E = 14 MeV and T = 13 MeV. 

What can one say about possible resonances in other A = 3 nuclei? The effect of the 
Coulomb force in the mirror three-proton system is trivial, pushing the | resonance toward 
higher energy and increasing its width accordingly. The resonance parameters are E = 15 
MeV and T = 14 MeV. 

The situation is more complex in the cases of 3 H and 3 He, because these nuclei have 
two-body channels below the three- nucleon thresholds. Recently an experiment, studying 
the H( 6 He,a) reaction, found indications of an excited state of the triton at E* = 7.0 ± 0.3 
MeV excitation energy with T = 0.6 ± 0.3 MeV width ||. This state would be situated 
between the d + n and n + n + p thresholds. It was suggested in || and Jl7| that the 
dominant configuration in this state is the (nn)p one, i.e., the resonance originates from the 
closed three-body channel. 

Unfortunately our current N — N interaction is not appropriate for a calculation which 
contains all the relevant angular momentum configurations for the triton. This is because 
the Minnesota force reproduces the 3 S\ n + p phase shift and the deuteron binding energy 
without the tensor force. However, we can estimate the role of the n+n+p and d+n channels 
in a restricted model, which contains 3 Di states just between the d and n but not inside the 
deuteron. Th Minnesota force is shown |1| to give a good overall description for the d + p, 
and supposedly for the d+n, scattering in such a model. It is important to note that, if there 
is a resonance in the (nn)p configuration, it must appear in a one-channel d + n model, too. 
The reason is that any J n state that can be formed in the (nn)p configuration can be built 
up from d + n one, too, and this d + n channel is nonorthogonal to the (nn)p configuration. 
Thus, as the CSM is not well suited for finding poles at energies with negative real and 
imaginary parts (the E* = 7.0 MeV energy would mean such a position in the n + n+p 
channel), we look for poles in the d + n scattering. We use the analytic continuation of the S 
matrix to complex energies [|14||, which is feasible for a two-body case and equivalent with the 
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CSM. A i + pole is found at E = 1.4 MeV d+n center-of-mass energy with a width of 9 MeV. 
The n + n + p channel can be taken into account approximately, by using a few discretized 
continuum state of (nn) with positive energies and square-integrable wave functions. The 
inclusion of this n + n + p channel does not change the pole position significantly. The 
dominant configuration of the wave function is d + n, in contrast to the case of Jl7| . 



The quality of our 3 Si and 3 D\ N — N forces makes it impossible too reach a firm 
conclusion concerning the existence of a narrow excited state of the triton. Our main result 
is that the picture of || and [[H]], namely that this state originates from the (nn)p channel, 
is questionable. As such a state would be nonorthogonal to the d + n channel, this latter 
configuration must be important. 

Finally, we note that the use of realistic interactions together with a full model space 
(including I > 3) may change the parameters of the | + resonance. This is why we have 
made no attempt to further optimize our basis (expansion length, Gaussian widths, etc.) 
and find more precise resonance parameters. Our aim was to show the existence of such a 
state. 



IV. CONCLUSION 

In summary, we have searched for three-neutron resonances using the complex scaling 
method. We have found a J 77 = | resonant state at E = 14 MeV energy with T = 13 MeV 
width. All other partial waves up to J = | were found to be nonresonant. The parameters 
of the mirror | three-proton resonance are E = 15 MeV and T = 14 MeV. Our results are 
inconclusive regarding the narrow excited state of the triton, but show the importance of 
the d + n channel. 
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FIGURES 



FIG. 1. S-wave (a) and P- and D-wave (b) p + p phase shifts in the center-of-mass frame. 
Experimental data are taken from [15]. 



FIG. 2. Energy eigenvalues of the complex scaled Hamiltonian of the (a) I , and (b) | 
three-neutron states; E = Re(e), T = 2 Im(e). The dots are the points of the rotated discretized 
continua, while the circle is a three-neutron resonance. The rotation angle is 0.4 rad. 

FIG. 3. The 9 trajectory of the | + three-neutron resonance with a two-component wave func- 
tion. The 9 values are indicated. 
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